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In low density regime , the fluorescence of Frenkel exitons in crystal slab can be studied 
without the aid of rotating wave and Marckoffian approximation. The equations for the 
, case of double and triple lattice-layers are now solved exactly to give the eigen decay rates, 

' frequency shifts and the statistical properties of the fields. 
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I. INTRODUCTION 

^ I The fluorescence of excitons in a quantum well or crystal slab is of charater of collective radiation, since 

'n ■ exciton, as exited state of the whole quantum well (crystal slab) has a collective transition dipole moment. 

I However not all these collective radiation are superradiance. Actually, the exciton emission has many eigen 

modes, some of them are superradiant modes and some of them are subradiant modes. 

Exciton have important application in photonic devices because excitonc device may have small size, 
, low power dissipation, high speed and high efficiency. All of these are needed by integrated photo-electric 

' circuits. 

It is well known that the exciton in bulk crystal does not radiate, but forms polariton insteadl^'^1. This 
^ • shows that a general treatment of exciton radiation should take the reabsorption effect into account. 

I I There has been quite a lot of theoretical studies on fluorescence of excitons['^~^°l. In the case of Frenkel 

^ . excitons of low density, Knoester studied'^! the crossover from superradiant excitons to bulk polaritons 

' when the number of the lattice-layers in the crystal increases without bound. He gives the correct form of 

^ , the eigen equation for the frequency shift and decay rate in the single lattice layer case. We have pointed 

O^' outi"! that it should count in the two-photon coupling term properly in order to get this correct form 

of eigen equation. Knoester proposed'^] that F^k' which describe the coupling between the excitons of 
wave vector k and k' by exchange of photons is strongly peaked around k = k' and one may keep only 
, ^ the diagonal elements to a good approximation. There are some ambiguities in this proposition, since the 

^ . set of the values for k exists different selections. Besides, when one derives the whole set of eigen decay 

rates, some of them are small (subradiant modes), the contribution from off diagonal elements could be 
important. 

We have studied the single lattice layer case in some detail^^. However in the case of Frenkel exciton 
fluorescence, neither single layer case nor the very thick case is important in practice. In this paper we 
will study the fluorescence of Frenkel exciton in thin crystal fllm of double and triple lattice-layers. The 
Heisenberg equations without rotating wave approximation are solved in the low density regime without 
Mackofflan approximation. The two-photon coupling term in the interaction Hamiltonian is included 
properly. All eigen decay rates, frequency shifts as well as the evolution of fields in terms of their initial 
values are obtained consequently. We note the non-diagonal elements of Fkk' is essential to derivation of 
these results. Our approach can be readily generalized to the case of more lattice layers. 
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The exciton-phonon interaction is not taken into account in this investigation. ° We shall give a brief 
review of general formulation in section 2. Section 3 and 4 are devoted to the cases of double and triple 
lattice-layers respectively. A brief conclusion is given in section 5. 



II. BRIEF REVIEW OF GENERAL FORMULATION 

We first write down the general formulation for the crystal slab of N lattice-layers. The crystal is 

assumed to have simple cubic structure. When the two-photon coupling term -^^A is taken into account, 
the interaction Hamiltonian between exciton and photon for the low density excitons are described as'^^l 

Hint = nJ2G{q)0{k + q)[Bk{t) + Btk{t)][a,{t) + algit)] + 

q,k 

% 1g(«)G(5')0(5' -fc)0(fc + a)KW +ai,W][a,,(i) +al^,(t)] (1) 

where q and k are wave vectors for photon and exciton respectively, they are in the z-direction, perpendic- 
ular to the crystal slab, Bk{t) and aq{t) are exciton and photon annihilation operators respectively, Bk{t) 
and B'^{t) are assumed to satisfy the boson commutation relations 

[Bk{t), B+{t)]=6kk', (2) 
and O is the frequency of the isolated lattice atom, G{q) is the coupling constant, 

where V = AL is the normalization volume for the photon, Nt is the total number of lattice sites in the 
crystal s\a.h{NT = NNl, Nl is the number of lattice sites in each layer), d denotes the transition dipole 
moments of a single lattice atom supposed for simplicity perpendicular to the ^r-axis. 
In Ref . [5] it is said that k takes the discrete value 

7 27rTO , 

with m = 0, 1, • • • A'' — 1. But in the eq(l) of Ref.[5] as well as in our Hint, ^'s are assumed symmetric with 
respect to the zero, hence we should take m= —^{N— 1), - ■ ■ ,^{N — 1) instead'^^l. The values of q is as 
usual: 



27r? 

9=-^, j=0,±l,±2,.--,±oo. 
0{k — q) is the wave-vector matching factor, now takes the forml^^l 

0(.-,) = i "i:" e'C-*^! -"^'*-'""" , (5) 
A f-^ A s\a^{k-q)a 

I is the index for the layers. In case k = q, 0{k — q) will equal to one, when k ^ q and A' is sufficient large, 
0{k — q) will be small. 

We notice that no rotating wave approximation is made in eq.(l), and the two terms on the right-hand 
side correspond to single-photon coupling and two-photon coiipling respectively. 

From Hint and the commutation relation, one may immediately write down the Heisenberg equations 
for Bk{t), dq{t) and their hermitian conjugates. These equations are linear equations so that they can be 
solved exactly. Carrying out the half-side Fourier transformation 



2 



Bk{^)= / Bfe(t)e*'^*di, B+^iuj)= / B_k{t)e'^*dt, etc. 

Jo Jo 

as did in Ref.[5] and then eliminating the photon operators, we get 

o, .2 

k' 

=i[{uj+ n)Bk{o) -{^- ^)B+i^m - 2i^Y.^kk'{uj)[Bk'{o) - B+^,m 



k' 



where -8^/(0) means -Bfe/(i = 0), etc, and 



= E - q)0{k' - q) 

g 

Nailp f, 0{k - q)0{q - k') 



I 



(6) 



A 2 / ^« 2 27^ ' (^) 

in which 

Note that our Fkk' (t^) and 0{k — q) arc somewhat difFcrcnt from those defined in Re[5]. For simplicity we 
have negelected the term of static dipole-dipole interactions 

We should mention that even disregarding the terms proportional to 0^(0) and a_g(0), eq.(7) is still 
different from the result of Ref.[5] as mentioned in Ref.[ll]. 



III. THE CASE OF DOUBLE LATTICE-LAYERS 



Now we consider the special case N = 2. The wave- vector matching factor for A'' = 2 becomes 

0{k - q) = cos^{k - q)a (10) 

where k takes the values Wc shall use F++{uj), F^-{uj) to represent jl^+ji. (w), F-^-^{w) etc. 

and evaluate the integrals in eq.(8) by contour integration. The results are 



F++(a;) = F__(a;) = -i^, = = -i^e^ (11) 



where 



We note that the nondiagonal elements {F^ and -F-+) is of the same orderly as the diagonal elements 

{F++ and F__). 
The coupled equation(7) now becomes 

{J^ + iiori - n'^)[B+{uj) - B+{uj)]+iijorie^[B-{uj) - B+(w)] = Ao{w), (13a) 
iujr]e^[B+{ui) - Bt{uj)] + (w^ + iwr? - ^'^)[B_{uj) - B+{ui)] = -A^{-oj) (13b) 

where 
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Ao{u) = i[{LU + n + iv)B+{0) -{Lu-n + i-q)I3+{Q)] - r/e^ [B-{0) - B+{0)] 

A+{-Lu) = i[{uj -Q + iv)BX{0) - (w + O + ir?)B-(0)] - rye"^ [B^{0) - B+{0)] 

-2V2.^i:G(,)cos(f + f )[^ + 4#r]- (l^d) 
4 2 LJ — \q\c uj + \q\c 

Eqs.(13) are easily solved to get: 

- M = {u^' + ^u.V-n^)Ao{u;)+^u;ve^AU-u>) ^^^^^ 

{lo'^ + iujrj — + iujrje <= ){uj'^ + iajr] — ^l"^ — iurje ) 

B_M - Bt{^) = ^u.^e^Mu) + {^+^u^-n^)At{-u:) ^^^^^ 

(w^ + iLori — fi^ + iuj-qe <= ){uj'^ + iur] — fi^ — iwrye <: ) 



The roots of 



and 



u)^ + iuT] — + iwrye'"" = (15) 



+ iurj — Q^ — iurje*"" = (16) 

will determince the eigen decay rates and corresponding frequency shifts. These roots can not have positive 

imaginary part, namely the poles of 13^{uj) — B^{oj) will not be localized in the upper half w-plane, since 
two necessary conditions can be deduced for eq.(15) to have root of positive imaginary part: 

r]>2Q. and — > 27r 
c 

and both of these conditions are untenable. Similarly eq.(16) also can not have root of positive imaginary 
part. These results mean the basic physics laws will not be violated as in the case of monolayer I-^-^l. 
We have derived four physical roots of eqs.(15) and (16) as: 

u)i = rii — iFi, CO3 = —ill — *ri 

IV2 = ^2 — i^2, W4 = —Q2 — (17) 



in which 



n2 = n{i-^£) , r,--,:-^. (is) 

In the following we will omit the terms proportional to 0^(0), a^(0) as did in Ref[5], since here we just 
study the fluorescence of excitons. Then the electric field is derived as foUowsI^^i 



1 n^a^ 
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oo+ie 



E{z,t)^--I du)E{z,Lj)e-^\ (19a) 
27r 



E{z,uj) = ^^*^[a,(a;) - a^,(a;)]e^«^ (19b) 
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The summation of q in eq.(19b) can be tranformed to integration and carried out by contour integration. 
In the positive z region outside the crystal slab, we get for the double lattice-layer case 

, .2Trnd, Lja{iv + n){B+{0)+B_{Oj) + {iv-n){B+{0) + B+{0)) 

bj[Z,LO) = I T=[COS — j;^ 

ac\/ A 2c ui"^ + iujr] — Q!^ + iuir}e~ 

. (g; + ^^)(^+(0) - ^-(0)) -{u^- O)(g|(0) - B+{Q)) ^ 
+ lujr] — \V — iu)r]e <= 

where A is the area of each layer, it is also the cross area of the normalization volume for the photon as 
mentioned above. 

The electric field E{z,t) in this region is calculated by eq.(19a), with the results 

E{z,t) = Q, for z-ct>0 (21a) 
E(z,t)=s{z,t) + h.c., for z - ct < (21b) 

where 



*^ = -cy^^' + n;-'n;^''' — 2^ — 

[(^+(0) + B_(0)) + ^li-ll-£i(B+(0) + Sl(0))]e-"^(*-^)-r^(*-^) 

+ -c\l^^' + n-r'n-2^'''' 

[(S+(0) - S_(0)) - ^i^§^(S|(0) - 5l(0))]e-^^(*-f ). (22) 

We note that this solution is free from Marckoffian approximation and also free from antirotating wave 
interaction. 

The electric field in the 2; < region can be derived similarly, with the resultant waves propagating in 

backward z direction as expected. 

There are two eigen decay rates appeared in the E{z,t): Ti and The corresponding eigen modes 
are linear combination of the two modes of m = ^ and m = — 5. As can be seen from eq.(20), these two 

eigen modes, which will be called as supcrradiant mode and subradiant mode , correspond to the operators 
-^{B+{0) +i3_(0)) and ■^(i3+(0) — i?_(0)) respectively. Hence, They correspond to modes of fc = and 
k = 1 with the operators 

Mt) = \ E ^'"""'Mt) = 7s[^+(^) + B-{t)], (23a) 



Bi{t) = \ E e^(^-5)'«B,(i) = ^[S+(t)-S_(t)]. (23b) 

Evidently, the dipoles of the two layers have the same phase for the former and have opposite phase for the 
latter. We note that the decay rate r2 of subradiant mode here is still as large as ^ times the decay rate 
of a single atom(molecular), because the atoms in each layer are still cooperated. The decay rate of = 
mode is twice of that of monolayer, which is just the character of superfluorescence. As can be seen from 
eq.(22) that even for the superradiant mode in which the emission is totally collective, the emitted light 
still may have different statistics and coherent properties according to the initial exciton state (also see 
the discussion in Rcf[ll]). For example, the coherent part of the electric field < E{z,t) > will be nonzero 
when the initial state of the exciton is a coherent state. But when the initial density matrix of the exciton 
is diagonal in Fock representation (including number state, chaotic state), the coherent part of < E{z, t) > 
will be zero. 

Up to the first order of ^ and ^, the E{z,t) is expressed by the superradiant mode operator Bq{0) 
and the subradiant mode operator Bi(0) as follows: 
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+ /%^(^)^^We--^^*-^-^(*-^+/..c. (24) 

for z > and t - ^ > 0. Similar results for ^ < 0, i + ^ > 0. E{z, t) is equal to zero if > , t - < 0) 
ov {z <Q,t+ I < 0) 

Since we have seen irregular behavior in the usual intensity operator for the solution of single layer 
casein], here only the energy flux operator is given instead. The energy flux is usually defined by 

S{z, t)^-^ : t{z, t) X B(2, t) : . (25) 
47r 

It is readily to show that S is always directed outward from the crystal film. So we rewrite S as nS 
which n is unit vector directing outer space from lattice-layers. Namely, it is in positive z direction in the 
.2 > region and in negative z direction in the z < region as required. 

So we may obtain S{z,t) from eq.(24). After neglecting oscillating terms and higher order terms of ^ 
and ^ (only keep first order terms), we have: 

S{z,t) = ^[B+(0)Bo(0) + ^BiiO) - i^So+^(0)]e2''(f-*) 
+ !Z^4+(0)Bi(0)e2'''(f-*) 

+ ^^^^[B+mm + Bt moio) 



^(Bo+(0)Bi(0) + Bo(0)Bi(0) 
- B^(0)Bo(0) - 4+(0)Bo+(0))]e(''+''')(t-*), (26a) 



with 



V' = vi^. (26b) 

We see from eq.(26) that the energy fiux decays in three different rate. The first term which is con- 
tributed by the cxciton of the short lifetime^ palys a important part at the bcgining time. The second term 
contributed by the exciton of the long lifetime becomes dominat at late time. The third term will exhibit 
itself in the intermediate stage. 



IV. THE CASE OF TRIPLE LATTICE-LAYERS 

The cases of odd N and even N have a qualitative difference in the m- value series —^{N—1),---^{N—1) 
for cq.(4). In the former case, m contains zero, while in the latter, not. = 3 is the simplest case of odd 
N, apart from the trival case iV = 1, which has no nondiagonal terms Fmm' (here and in the following 

27rm 
Na 



we use Fmm' to denote Fkk' according to the relation k = 2^). Thus we will study it as an example. For 



3 

leading to the matrix F (with elements Fmm', m,m' = 1, 0, —1) as 



0{k -q) = h2 cos(fc - q)a + 1] (27) 



F{uj) = -^-^ — I -x^+x 2x2 + 4x + 3 -x^ + x \ = -^r]D{uj), (28) 
2x^-2x -x^+x 




where x = e'"° = e*''. To the second order of 6, 
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A5 + 3i(52 ,5 + |i(52 -25 - SiJ^ 
D{uj) =1 5+ 9i - 8^ - 6i(52 5 + fiJ^ | . (29) 

-25 - 3i^2 ^ + |i^2 4^ + 3^^2 

We see that the nondiagonal elements are of the same order of Fn and F_i _i, so that they can not be 
neglected in the equations for Bi — B'ti and B-\ — B^ {Bm also means Bk for k = ^jf^)- 
The couple equations now take the form 

m' 

= i[{uj + n)B^{0) -iuj- il)B+^m + ^r?^Z?wM[^„'(0) - B+„,(0)] (30a) 



with 



m,m' = -1,0,1. (30b) 



In eqs.(30) the terms proportional to 0,(0) and alg(O) are neglected. 

The eigen decay rates and corresponding frequency shifts are determined by the roots of the following 
equation 

(^2 - f22)3 ^ 1 ^ 2A){w^ - fi2)2 + LW{2AB + A^-E^- 2C^){lo^ - fi^) 

+ ^ifuj'^[{A^ - E^)B + 2C^{E - A)] = 0, (31) 

where A, B, C and E are matrix elements of D{lo), defined as follows: 

( A{i^) C(w) E{u) 
D{u}) = C{oj) B{uj) C{uj) I . (32) 
\E{co) C{oj) aIoj) 

We get the six roots of eq.(31) as follows: 

wi = Oi — iTi, co[ = — Oi — iVi, 

uiQ = ^ — *r, ~ ^ 'To, (33) 

ui-i = Q-i — iT-i, u>'_i = —Q-i — iV-i. 



with 



fii^fi(l-l^) , r^ = j^n"-^, (34a) 

0_i = (l-^) , r_i=r?^. (34c) 

All the roots are in the lower half plan of complex oj as they should be. 

The direct way to solve for Bm{tjj) — B^^{co) from eq.(30) is to diagonize the matrix D{u)) defined by 
eq.(28). Up to second order of 5, wc get the transformation matrix T{co) which satisfies 

T{u))D{w)T{u}) = I £>o(w) I (35a) 



1 Q2a2 



as 



r(c.) 



/ M_ -V2CM M 

' V2 B-A-E -J2 

M ^ I , (35b) 

\ V2 ^ ~V2 
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where 

M = ^ . (35c) 

The result for -8^(0;) — 5+^(0;) is expressed then by 

^rm',m ^2_o2+ 1 , [('^ + ^ + ^^^™')/3i^^('^,0) - (a; - n + ^r?i?„0C' ('^,0) (36a) 
in which 

/3(,i) (w, 0) = J 5^ T„™' (a;)B„. (0), 0^^^ {u, 0) = ^ ^ T^™' (a;)Bi„, (0) (36b) 

Substituting eqs.(36) into eqs.(19) and carrying out the integrations, it finally yields 
aq{u) - atgico) = -2-^^^G(g) ^ [2cos(^ - qa) + l]x 

m,m' 



a;2 _ 02 + i?7w£>„ 



m + ^Wr:,! 0)-in- iv)p'^> {w, 0)], (37a) 



m.m' 



^2 _ f]2 ^ IrquDrt 



m + ^WJ (a;, 0) - (n - a;)/?;^^/ (c., 0)] (37b) 



and 



Eiz, t) = V^e(t - ^) 5: d„e-^"-(*-f )-r'"(*-f ) + /..c, (37c) 



for ^ > 0, t - - > 0, where 



1 v^rn ,27rTO' WtoOx , X 

= ?^Z^t2cOS(— ) + Ijrmm'(Wm) 

[(fi + a;„)^^) (a;„, 0) - (fi - a;„)/3^) (a;„, 0)]. (37d) 
To the leading term, ctm are given by 

ai=i^[Bi(0)+B_i(0)], 

do = 4(0), (38) 
d_i = ^[Bi(0)-^_i(0)], 

We see that m = ±1 modes are not of eigen decay rates. On the contrary, the eigen modes are nearly 
maximum mix of these two modes. 

In terms of the creation operator for an excitation in the /th layerl^'^^1, we have 
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Thus we have, denoting Bk by as before, the three eigen modes. They are approximated as following 
;^[^„=i(0) +B„=_i(0)] = -^[-Bi=i(0) + 2Si=o(0) - B+^_,], (40a) 

Bm=o = -^[Bi=m + Bi=o{0) + -B(=-i(0)], (40b) 
v3 

-^[Bm=m - Bm=-im = -^[Bi^o) - Bi=-im, (40c) 

while 

Bm=±m = ^[{-l T ^i)Bi=m + Bi^oio) + {-I ± y^m^-im- im 

We see the superradiant mode(m = 0) has a decay rate 2ro = 3r] which is triplet of that for mono- 
layer, showing the emission is totally cooperative. However, as in the double layer case, the statistical 
properties of the light of this mode still may have different varieties which depend on the initial state of 
the excitons. For the z > a region, we may also give the energy flux of the case of triple lattice-layers 
according to eq.(25) and eq.(37) as following: 

{S{z,t)) = {S,{z,t)) + {S2{z,t)), (41) 

where {Si{z,t)) is the main part, it is expressed by 

{Si{z,t)) = ^{9(B+(0)5o(0))e-3''(*-f) 



9c2 



-(B+(0)B+(0))e-^ 



{S2{z,t)) = -t-^^[{B+{0)Bo{0) - B+{0)B+{0)) 



+ ^((Sl(0)B_(0))e-«''(*-t)} (42) 
with B±{0) = -^(51(0) ± B_i(0)) and {S2{z,t)) may be approximated by 

i3V3{Bl{0)Bo{0) - Bo+(0)B_(0))]e-i''(*-^). (43) 
It only appears when the initial exciton density matrix in Fock representation has non-diagonal elements. 



V. BRIEF SUMMARY 

1. KnoesterPl claimed that Fkk'iuj) is diagonal to a good approximation. But we show explicitly that 
this is not generally true. 

2. In low-density case, even the emission is superradiant in nature, the light still may have different 
coherent statistical properties, depending on the initial state of exciton. 
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